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ABSTRACT 

We carry out a comprehensive investigation comparing the three-dimensional magnetic field restructuring, flare 
energy release, and the helioseismic response of two homologous flares, the 2011 September 6 X2.1 (FL1) and 
September 7 X1.8 (FL2) flares in NOAA AR 11283. In our analysis, (1) a twisted flux rope (FR) collapses onto 
the surface at a speed of 1.5 km s _1 after a partial eruption in FL1. The FR then gradually grows to reach a higher 
altitude and collapses again at 3 km s _1 after a fuller eruption in FL2. Also, FL2 shows a larger decrease of the 
flux-weighted centroid separation of opposite magnetic polarities and a greater change of the horizontal field on 
the surface. These imply a more violent coronal implosion with corresponding more intense surface signatures in 
FL2. (2) The FR is inclined northward and together with the ambient fields, it undergoes a southward turning after 
both events. This agrees with the asymmetric decay of the penumbra observed in the peripheral regions. (3) The 
amounts of free magnetic energy and nonthermal electron energy released during FL1 are comparable to those of 
FL2 within the uncertainties of the measurements. (4) No sunquake was detected in FL1 ; in contrast, FL2 produced 
two seismic emission sources SI and S2 both lying in the penumbral regions. Interestingly, SI and S2 are connected 
by magnetic loops, and the stronger source S2 has a weaker vertical magnetic field. We discuss these results in 
relation to the implosion process in the low corona and the sunquake generation. 
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1. INTRODUCTION 

There is often a coincidence among flare emissions, launch of 
the coronal mass ejections (CMEs), and the impulsive seismic 

waves (“sunquakes”) embarking on and beneath the Sun’s sur- 

face, and they are usually considered to be related to each other. 
However, the transport of energy and momentum into the inte- 
rior from the solar atmosphere is not yet fully understood. Since 

solar eruptions are generally believed to be a result of magnetic 
reconnection in the corona (e.g., Forbes 2013; Su et al. 2013), 
changes in the magnetic field due to coronal field restructuring 
must be considered not only in studying flares/CMEs but also 
in understanding the coupling between the corona and the solar 
interior (Hudson 2011). Efforts in this regard were originally ini- 
tiated from the theoretical viewpoint. In particular, the concept 

of “implosions” in coronal transients (see outward explosions 

of CMEs) was put forward to predict an inward contraction of 
coronal field that would occur simultaneously with the magnetic 
energy release (Hudson 2000). Clear evidence of implosions 
during the flare impulsive phase has recently been gathered 
from time sequences of high-resolution coronal images (e.g., 
Liu et al. 2009, 2012c; Gosain 2012; Simoes et al. 2013). The 
implosion picture may be considered in line with the reconnec- 
tion model of two current-carrying flux loops, which produces 
two new current-carrying flux systems with one moving closer 
to the surface and another moving upward as part of CMEs 

(Melrose 1997). This kind of loop-loop interaction is depicted 
in the tether-cutting reconnection model for sigmoidal active 


regions (ARs; e.g., Moore et al. 2001; Liu et al. 2007, 2010, 
2013a). The downward drop of coronal currents after eruptions 
has also been indicated by comparing the pre- and postflare coro- 
nal field models reconstructed using the nonlinear force-free 
field (NLFFF) extrapolation technique (Liu et al. 2012a; Sun 
et al. 2012). 

A direct consequence of coronal implosions would be a more 
horizontal (i.e., more inclined to the surface) configuration of 
the photospheric magnetic field. This is supported by the finding 
that the transverse magnetic field around the magnetic polarity 
inversion line (PIL) at the center of the flaring region often 
exhibits a rapid and persistent enhancement immediately after 
flares/CMEs, as revealed by vector magnetograms of either 
ground-based or space-borne instruments (e.g., Wang & Liu 
2010; Liu et al. 2012a; Wang et al. 2012b, 2012c; Sun etal. 2012; 
Petrie 2013). Intuitively, such an inward collapse of the central 
magnetic field may be accompanied by an upward turning of 
the peripheral field in the outer areas of flaring regions. At these 
locations the photospheric transverse field would subsequently 
decrease, as sometimes reported (Li et al. 2011; Sun et al. 
2012). Since sunspot penumbrae in white light (WL) indicate 
an inclined magnetic field structure, it is remarkable to note 
that the above surmised collapse-and-turning scenario of flaring 
magnetic field (as simply but attractively illustrated in Liu 
et al. 2005) is well corroborated by the observed flare-related 
strengthening (decay) of central (outer) sunspot penumbra of 8 
spots (Wang et al. 2004, 2012a, 2013; Liu et al. 2005; Deng 
et al. 2005, 2011; Chen et al. 2007). 
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Importantly, the change of photospheric magnetic field at the 
flare core region to a more horizontal state implies an upward 
Lorentz-force change exerting on the outer atmosphere, which is 
balanced by an equal and opposite Lorentz-force change acting 
on the photosphere and below (Hudson et al. 2008; Fisher et al. 
2012). The latter is described as a “back reaction” caused by the 
coronal field reconfiguration (Hudson et al. 2008). It is asserted 
that the integration of the upward Lorentz-force change over its 
change period should produce an upward impulse driving the 
erupting CMEs; the conservation of momentum also signifies 
an equal, downward-moving impulse, which is applied into the 
solar interior and could supply enough energy for seismic waves. 
However, the flare-CME-sunquake relation as implied above is 
not yet fully understood. More particular conditions may have 
to be fulfilled to generate a CME or a sunquake since not every 
flare is associated with them. 

Sunquakes are produced in a response of the low solar 
atmosphere to the localized hydrodynamic impact of flares/ 
CMEs. As a result, acoustic waves are generated, and they 
travel into the interior and refract back to the solar surface 
(for reviews, see Donea 2011 and Kosovichev 2014). Although 
circular ripples emanating from a flare site were discovered 
more than a decade ago (Kosovichev & Zharkova 1998), the 
physical origin of this wave phenomenon is still unknown. 
The main driving mechanisms of sunquakes proposed thus far 
include precipitation of high-energy electrons or protons (e.g., 
Kosovichev 2007; Zharkova & Zharkov 2007), back-warming 
radiation heating related to WL emission (e.g., Donea et al. 
2006; Donea 2011), and changes of magnetic field (e.g., Hudson 
et al. 2008; Martmez-Oliveros & Donea 2009; Zharkov et al. 
2011; Fisher et al. 2012; Alvarado-Gomez et al. 2012). As the 
production of seismic activities could involve multiple layers 
of the solar atmosphere, study of sunquakes will advance the 
current understanding of solar flare physics. 

We note that the coronal implosion scenario includes multiple 
components for back reaction, and the role of each component 
has been studied separately. In this paper, we carry out a com- 
prehensive investigation of the homologous 2011 September 6 
X2. 1 and September 7X1.8 flares, including the coronal field re- 
structuring using the NLFFF modeling, photospheric magnetic 
field and sunspot structure changes, energy release during flares, 
and the helioseismic response to the flare impact. We consider it 
instructive to inspect the differences between these two events, 
because such a comparative study could be important to better 
explore and understand the role and relationship of different 
phenomena associated with flares. The plan of this paper is as 
follows. In Section 2, we first provide an events overview and 
then introduce the observations and the data processing proce- 
dure. In Section 3, we describe the observational and model 
results and point out their implications. Major findings are sum- 
marized and discussed in Section 4. 
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to the flares. They found a significant increase of the horizon- 
tal field component that corresponds to a notable Lorentz-force 
perturbation. These changes are considered rapid (in ~30 min- 
utes) compared to the long-term evolution of the AR field on 
a timescale of days. Jiang et al. (2014) examined the buildup 
process and instability condition of the AR sigmoid before the 
X2.1 flare using a NLFFF extrapolation method. A realistic ini- 
tiation of the eruption was then achieved by Jiang et al. (2013) 
based on an MHD model. Feng et al. (2013) analyzed the en- 
ergy budget of the first event, and concluded that this X2.1 flare 
and the associated CME had similar total energies. The sunspot 
rotation before the onset of the X2. 1 flare was also observed and 
linked to the event onset (Ruan et al. 2014). In particular, the 
second flare (XI .8) produces a pronounced sunquake as reported 
by Zharkov et al. (2013); however, no clear seismic signatures 
were detected during the earlier X2.1 flare, and these homol- 
ogous flaring events seem to have similar emission patterns in 
WL and X-rays (Xu et al. 2014). 

In this analysis, full-disk photospheric vector magnetograms 
are supplied by the Helioseismic and Magnetic Imager (HMI; 
Schou et al. 2012) on board the Solar Dynamics Observatory 
(SDO\ Pesnell et al. 2012). The vector field data are derived 
using the Very Fast Inversion of the Stokes Vector algorithm 
(Borrero et al. 2011), and the minimum energy method (Metcalf 
1994; Leka et al. 2009; Metcalf et al. 2006) is used to resolve the 
180° azimuthal ambiguity. We acquired vector magnetograms 
of AR 11283 that were prepared by the HMI team for special 
campaigns for selected ARs. 9 The images are remapped using 
Lambert (cylindrical equal area) projection centered on the 
midpoint of the AR, which is tracked at the Carrington rotation 
rate. The observed fields are also transformed to heliographic 
coordinates (Gary & Hagyard 1990) to remove the projection 
effect. The vector data analyzed have a pixel scale of ~0'. / 5 
and a cadence of 12 minutes, spanning from 20:22 UT on 
2011 September 6 to 00:46 UT on September 8 throughout 
both the X2.1 and X1.8 flares. The former (latter) started 
at 22:12UT (22:32 UT), peaked at 22:20 UT (22:38 UT), and 
ended at 22:24 UT (22:44 UT) on September 6 (7) in GOES 
1-8 A flux. 

Modeling of the coronal field was carried out using NLFFF 
extrapolations. The net force and torque in the observed photo- 
spheric field were first minimized by a preprocessing procedure 
in order to obtain the chromosphere-like data that meet the force- 
free condition (Wiegelmann et al. 2006). The “weighted opti- 
mization” method (Wheatland et al. 2000; Wiegelmann 2004) 
was then applied to the preprocessed photospheric boundary to 
extrapolate the magnetic field toward coronal heights. This ap- 
proach minimizes an integrated joint measure L over the volume 
V, which comprises the normalized Lorentz force, the field di- 
vergence, and treatment of the measurement errors (especially 
for the transverse field B T ; Wiegelmann & Inhester 2010): 


2. EVENTS OVERVIEW AND DATA PROCESSING 

The source NOAA AR 11283 emerged on the east limb 
on 2011 August 30. No major flares occurred until Septem- 
ber 6, when the region appears in the /3y / fid configuration 
containing one 8 spot (see Figure 1). The September 6 X2.1 
and September 7 XL 8 flares are the only X-class flares that 
occurred in this AR, and they make up a pair of homologous 
flares/CMEs due to similar locations and emission signatures 
in multiple wavelengths. Several studies of these flares have 
been conducted. Wang et al. (2012c) and Petrie (2013) studied 
the evolution of the photospheric vector magnetic field related 


-L 


|(V x B) x Bl 2 , , 

Wf- y L + Wrf |V -B| 2 d 3 V 


B 2 


+ V 


l (B - 


B obs ) w (B - B obs ) d 2 S, 


( 1 ) 


where B is the magnetic field vector, w f and Wd are weighting 
functions, v is a Lagrangian multiplier that controls deviations 
between the model and observed (B 0 b s ) fields, and W is a space- 
dependent diagonal matrix incorporating errors. For optimal 


9 http://jsoc.stanford.edu/jsocwiki/ReleaseNotes2 
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Figure 1 . HMI continuum intensity images right before the X2.1 flare on 201 1 September 6 (a) and the XI. 8 flare on September 7 (b), showing the evolution of NOAA 
AR 11283 in one day. The solid and dotted contours represent positive and negative magnetic field, respectively, at levels of ±200 and ±1000 G. p and f indicate the 
preceding and following spots, respectively. The images have a field of view of 67" x 67" (same as Figures 2, 4(a)-(d), and 8(a)-(f)). 


results, we set v = 0.001, the vertical (z) component of W to 
1, and the horizontal components (x, y) of W to 
(Wiegelmann et al. 2012). Using 2x2 rebinned magnetograms 
covering the entire AR, the field extrapolation was performed 
within a box of 300 x 256 x 256 uniform grid points, 
which corresponds to about 219 x 187 x 187 Mm 3 . Note 
that the bottom boundary is assigned a uniform altitude of 
0.73 Mm (1 pixel) in this study as a result of the preprocessed 
bottom boundary. For the 2 hr period before and after both 
the X2.1 and XI. 8 flares, the calculation was conducted at the 
full cadence (12 minutes); between the flares (from 00:22 to 
20:22 UT on 2011 September 7), a lhr cadence was adopted. 
To quantitatively evaluate the model quality, we resorted to 
a current- weighted (CW) average of sine metric (CWsin0) 
where 0 is the angle between B and the electric current J, 
and a mean absolute fractional flux change metric (\fi\) = 
(|(V • B)|/(6|B|/Ax)) where Ax is the grid spacing (e.g., 
Wheatland et al. 2000; Schrijver et al. 2008). These metrics 
have a zero value given a perfectly force-free and divergence- 
free field. For the inner computation box (excluding the buffer 
boundary of 32 pixels in the lateral and top) of our 63 NLFFF 
models, (CWsin6>) = 0.26 ± 0.02 and (\ft\) x 10 4 = 5.4 ±0.9, 
which are typical for this algorithm and suggest moderately 
satisfactory solutions. 

Evolution of the AR field can be revealed by surface flows, 
which were traced with the differential affine velocity estimator 
for vector magnetograms (DAVE4VM; Schuck 2008). A win- 
dow size of 19 pixels for feature tracking was set according 
to former studies (e.g., Liu et al. 2013b). HMI also produces 
45 s cadence continuum intensity images, Dopplergrams, and 
line-of-sight (LOS) magnetograms. The WL emission of flares 
as well as sunspot structure variation were examined using in- 
tensity images. A high-resolution ((X'l) TiO (a proxy for con- 
tinuum at 7057 A) image taken by the New Solar Telescope 
(NST; Goode et al. 2010; Cao et al. 2010) at Big Bear Solar 
Observatory (BBSO) further helps to recognize the sunspot WL 
structure in detail. Based on the Dopplergrams, the acoustic 
sources of sunquakes were probed by applying the helioseis- 
mic holography technique (e.g., Donea et al. 1999; Lindsey & 
Braun 2000; Donea & Lindsey 2005), which reconstructs the 
acoustic emission by producing phase-coherent snapshots of 
acoustic egression power. Time profiles of the integration of the 


source area on the snapshots then give the temporal evolution 
of sunquake sources. 

The hard X-ray (HXR) emission of the 2011 September 6 
X2.1 flare was registered by the Reuven Ramaty High Energy 
Solar Spectroscopic Imager ( RHESSI ; Lin et al. 2002), and the 
HXR images and spectra were investigated by Feng et al. (2013). 
The 2011 September 7 XL 8 flare occurred during RHESSI 
night but it was observed by the Fermi Gamma-Ray Burst 
Monitor (GBM; Meegan et al. 2009), with data analysis made 
possible by the Fermi Solar Flare Observations facility. 10 Fermi/ 
GBM produces HXR time profiles and spectra with no imaging 
capability. In addition, to identify coronal magnetic structures 
(e.g., filaments and hot coronal loops) related to eruptions, we 
used 304 A (He n; 0.05 MK) and 94 A (Fe xviii; 6.3 MK) images 
taken by the Atmospheric Imaging Assembly (AIA; Lemen et al. 
2012) on board SDO. 

3. RESULTS AND ANALYSIS 

In this section, we describe observational and model re- 
sults obtained using all data sets, and discuss their significance. 
Specifically, we first present the flare-related changes of mag- 
netic field on the surface, then proceed to the three-dimensional 
(3D) magnetic field restructuring (more dynamic detail can be 
seen in an accompanying animation in the online journal), and 
subsequently estimate the amount of the free magnetic energy 
released during the flares as well as that of nonthermal energy. 
Finally, we study seismic sources and their relation to flare 
emissions and magnetic field structure. 

3.1. Surface Magnetic Field Changes 

Both the 2011 September 6 X2.1 and September 7 XL 8 
flares are initiated around the highly sheared PIL of the AR 
(Figures 2(a) and (e)). In this area magnetic field concentrations 
with opposite polarity show a converging motion in the N-S 
direction and a diverging motion in the E-W direction, which are 
clearly tracked by DAVE4VM (Figures 2(c) and (g)). The main 
proceeding (p) and following (f) spots of the 8 configuration 
gradually separate from each other from September 6 to 7, 
which eventually “straightens” the PIL (also see Figure 1). 


10 http://hesperia.gsfc.nasa.gov/fermi_solar 
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Figure 2. Preflare HMI vertical magnetic field superimposed with arrows illustrating horizontal field vectors (panels (a) and (e)) and plasma flows (panels (c) and (g)). 
For clarity, arrows in negative and positive fields are coded in different colors. The flare-related changes of the horizontal field (panels (b) and (f)) and intensity (panels 
(d) and (h)) are presented. The intensity difference images are smoothed by a window of 275 x 275. The running difference images of intensity near the HXR peak are 
contoured at 3800 DN in panel (c) for the 2011 September 6 X2.1 flare (22:18:37 UT frame minus 22: 17:52 UT frame), and in panels (f) and (g) for the September 7 
X1.8 flare (22:36:37 UT frame minus 22:35:52 UT frame), showing flare impacts P1-P2 and S1-S4, respectively; see Section 3.4 for details. In the lower panels, the 
solid straight lines indicate the bottom side of the rotating vertical slice used in Figure 4. The dotted lines show the main segment of the PIL. 

(A color version of this figure is available in the online journal.) 


To better delineate the entire evolution and study the position 
change of the major photospheric magnetic field concentrations, 
we calculate the centroid (flux-weighted average) positions 
(x c , y c ) of both polarities as = J2 x i B r(i)/ J2 B r(i) and y c = 
J2 yi B Hi) / B r(i ) j where (xi , yt ) is the coordinate of pixels with 

an absolute vertical magnetic flux density \B r ^\ > Bq (Wang 
2006). For this AR, Bq is set to vary from 800 G to 1000 G for 
uncertainty estimation. The temporal evolution of the centroid 
separation between the main positive and negative fields is 
plotted in Figure 3(e), which demonstrates their separation 
motion at a speed of about 0.04 km s -1 at most times. It can 
also be obviously seen that the centroid separation decreases 
immediately after both flares (the vertical yellow shades) for 
a short time period, after which the long-term evolution trend 
is restored. Specifically, the separation of the two polarities is 
shortened by ~0.85 Mm in 1 hr (at ~0.26 km s -1 ) following the 
earlier X2.1 flare, and by ~1.4 Mm in 2.4 hr (~0.17 km s -1 ) 
following the later XI. 8 flare. These temporary but rapid 
approaching of the two components of the 8 spot are opposite 
to the overall shearing motion, implying the relaxation of the 
sheared magnetic field after the flares (Wang 2006) in accord 
with magnetic implosions (Hudson 2000). Here the decrease of 
the centroid separation distance associated with the XI. 8 flare 
is more significant (~1.6 times) compared to that associated 
with the X2.1 flare, implying a larger energy release in the 
former event. 

Rapid changes of the horizontal field B^ are also found in both 
flares. A stepwise increase of B^ (colored red in the difference 
B h maps in Figures 2(b) and (f)) is primarily observed within the 
elongated region lying along the central flaring PIL, as outlined 


using the box R1 (R2; also see Petrie 2013). There the mean 
#h abruptly increases by ~26% in the X2.1 and by ~38% 
in the XI. 8 flare (Figures 3(a) and (b)). Corresponding to the 
surface field change, the radial component of the Lorentz-force 
change acting on the surface and interior can be quantified as 
8F r = (1/87 r) / dA(8B^ — 8B^) (Fisher etal. 2012). Integrating 
over the region R1 (R2) yields a similar downward 8 F r peaking 
at 2.5 x 10 22 (2. 1 x 10 22 ) dyne during the X2.1 (XI. 8) flare. The 
increase of B^ around flaring PILs is a common feature of major 
flares (e.g., Wang & Liu 2010), and could be closely related to 
implosions of the coronal field (e.g., Liu et al. 2012a). 

More interestingly, the central areas of B^ enhancement are 
surrounded by the ring-like regions with decreasing B h (colored 
blue in Figures 2(b) and (f)). It is remarkable that the rings 
are not symmetric, but much more intense on the northern side 
in the fan-shaped sections Rid and R2d. The mean B h in Rid 
(R2d) tends to decrease for about 0.5 hr, with a sharpest drop 
by 12% (16%) in the X2.1 (X1.8) flare (Figures 3(c) and (d)). 

We note that (1) previous works are mainly focused on the 
increase of B^ in the flare core region, while the weakening 
of #h in the peripheral areas has not been studied in detail. 

(2) The distinct asymmetry of the decreasing B^ region is 
unlikely an artifact due to the projection effect. This is because 
the flaring AR is considered not far from the disk center (with an 
orientation cosine factor of 0.4 at the X2.1 flare and 0.54 at the 
XL 8 flare), and that both events show a similar change pattern. 
Instead, the asymmetric field change on the photosphere may 
connote an asymmetry in the 3D magnetic field structure and its 
evolution. This will be further discussed in the next section. 

(3) The ring-like change pattern of B h is very reminiscent 
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Figure 3. Temporal evolution of the mean horizontal field B h and inclination angle cp in the regions R1 (a), R2 (b), Rid (c), and R2d (d), as denoted in Figure 2, and 
that of the flux- weighted centroid separation of the two magnetic polarities of the vertical field (e). The green lines are the linear fits to the data. The yellow shades 
indicate the duration of the X2.1 and X1.8 flares in GOES 1-8 A flux. 

(A color version of this figure is available in the online journal.) 


of similar features in the flare-related darkening (decay) of 
the inner (outer) penumbra (Liu et al. 2005), which are also 
observed in the present events. In the difference images before 
and after the flares (Figures 2(d) and (h); see Figure 1), we 
see the central dark region (mainly penumbral strengthening) 
encompassed by asymmetric white areas (penumbral decay), 
which well corresponds to the regions of B^ increase and 
decrease, respectively. The magnetic inclination angle relative 
to the horizontal plane <p = tan -1 (|£ r \/Bh) also changes 
accordingly after both flares. It rapidly decreases by 14% (23%) 
in the R1 (R2) region, while increases gradually in the Rld/R2d 
regions by ~18% in 0.5 hr (Figures 3(a)-(d)). 

As will be shown below, we believe that the observed flare- 
related photospheric changes are an intrinsic manifestation of a 
reconfiguration of coronal fields. This view is consistent with 
the simulation of Rempel & Schlichenmaier (2011), in which 
sunspot penumbrae form or become destructed within ~0.5 hr 
after the boundary condition at the top is altered. 

3.2. 3D Magnetic Field Restructuring 

We first inspect the general coronal field structure right be- 
fore the X2.1 and XL 8 flares using our NLFFF extrapolations. 


Selected field lines within the main 8 spot region are plotted in 
Figures 4(a) and (c). These lines are colored according to the 
twist index T n = {\/An)otL = (iio/An)J r L/B r at their foot- 
points, where J r and B r are the vertical current and field on the 
photosphere, respectively, and L is the field line length (Inoue 
et al. 2011). It is clearly shown that sheared and twisted field 
lines lie along the PIL and are embedded in less twisted envelope 
field. The former core fields possess an overall sigmoid structure 
tracing out a flux rope (FR), which is manifested by a sigmoidal 
filament observed in AIA 304 A (Figures 4(b) and (d)). Com- 
pared with the X2.1 flare on September 6, more highly twisted 
fields up to 1.9 turns are present right before the XL 8 flare on 
September 7, together with an apparently thicker filament. This 
may imply that the XL 8 flare could be associated with a more 
well-formed FR (see further discussions below). Detailed topo- 
logical analysis of the FR as well as its instability condition are, 
however, out of the scope of this study (but see Jiang et al. 2014). 

In the following, we investigate the time sequence of static 
NLFFF models assuming that the coronal field evolution can 
be described by equilibria in successive steps (e.g., Regnier & 
Canfield 2006; Wu et al. 2009; Liu et al. 2012a; Sun et al. 
2012; Jiang et al. 2014; Jing et al. 2014). This approach allows 
us to gain an insight into the relaxation of the nonpotentiality 
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(e) IJ h l 1 1/09/06 22: 10:22 (f) IJ h l 1 1/09/06 22:34:22 (g) IJ h l 1 1/09/07 04:22:22 (h) IJ h l 1 1/09/07 10:22:22 



(m) a 11/09/06 22:10:22 


(n) a 11/09/06 22:34:22 


(o) a 11/09/07 22:22:22 


(p) a 11/09/07 22:46:22 



- 0.5 0.0 0.5 

Figure 4. Preflare HMI vertical field images superimposed with selected NLFFF lines (panels (a) and (c)), in comparison with cotemporal AIA 304 A images (panels (b) 
and (d)) overplotted with vertical field contours at levels of ±1000 and ±1500 G. The thick white slit in panels (a) and (c) indicates the bottom of a rotating vertical 
slice of 23 x 23 Mm 2 , over which the distributions of /h are plotted in panels (e)-(l) in logarithmic scale at selected instances. The white contours are at levels of 40% 
and 60% of the maximum current. The left bottom corner of the vertical slice corresponds to the southern end of the slit. The thick arrow in panel (e) indicates the FR 
orientation in the cross section. Panels (m)-(p) shows the distribution of the force-free parameter a in a larger slice (40 x 40 Mm 2 ) further extending northward, with 
the field of view of panels (e)-(l) denoted by the black box. The black arrows in panels (e)-(p) are the transverse field vectors in the vertical slices. See Section 3.2 for 
details. 


(An animation and a color version of this figure are available in the online journal.) 
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as well as the build-up process of magnetic energy between 
the eruptions. For our purpose of the comparative study of 
the 3D field restructuring associated with the X2.1 and XI. 8 
flares, we employ a vertical slice of 23 x 23 Mm 2 cutting 
through the middle of the FR and concentrate on examining 
the evolution of magnetic field and electric current in this cross 
section. The bottom side of this slice on the surface is overplotted 
on Figures 4(a) and (c) and also Figure 2 (bottom panels) as a 
thick straight line. The slice rotates counterclockwise at a speed 
of ~1?33 hr -1 so that it maintains an orientation perpendicular 
to the central PIL. In Figures 4(e)-(l) (and the accompanying 
animation), we depict the distribution of the horizontal current 
IJhl superimposed with the transverse magnetic field within 
this rotating vertical slice at selected instances, including those 
immediately before and after the flares. Some snapshots of the 
force-free parameter a = /roJ • B/|B | 2 in a larger slice are 
also presented in Figures 4(m)-(p). We describe these results in 
detail as follows. 

First, a spiral of field lines is evidently seen in the vertical 
cross section at most times, manifesting the poloidal flux of 
the FR. The system obviously carries free-magnetic energy as 
there exist significant field-aligned currents, with the maximum 
of Jk found to be generally cospatial with the axis of the FR 
(i.e., the center of the spiral; e.g., Canou & Amari 2010). The 
FR is also twisted and the maximum a near its axis reaches 
about 0.57 (0.48) Mm -1 right before the X2.1 (XI. 8) flare. Also 
noticeably, the FR is not symmetric in the vertical cross section 
especially before both flares, as it leans northward at about 66° 
relative to the surface. This is illustrated by the thick long arrow 
in Figure 4(e), which starts from the PIL and passes through the 
Jh maximum. The inclination of the FR toward north is mostly 
because that the magnetic flux (thus magnetic pressure) of the 
southern negative polarity is stronger than the northern positive 
polarity (Jiang et al. 2014). The overlying and neighboring field 
structure of the FR is delineated by the maps of the distribution 
of a in a wider field of view. 

Second, similar to the ascending motion of the FR before 
the X2.1 flare (Jiang et al. 2014), the FR gradually rises be- 
tween the occurrence times of the X2.1 flare on September 6 
and the XI. 8 flare on September 7 (Figures 4(f)— (j)). To aid in 
quantifying the FR evolution, we take the centroid of Jh as a 
proxy for the FR axis and follow its temporal evolution pro- 
jected onto the FR inclination direction. The statistical position 
error was estimated by varying the low threshold from 40% to 
60% of each maximum Jh (white contours in Figures 4(e)-(l)), 
as these contours cover a region similar to the observed fila- 
ment width. The result in Figure 5(a) shows that the altitude of 
the FR axis gradually increases from 2.7 Mm at 22:34 UT on 
September 6 to 7.1 Mm at 22:22 UT on September 7. During 
15-20 UT, the derived speed of this ascending motion has an 
average value of 0.09 km s -1 (see Figure 5(b)), which is com- 
parable to the slow rising speed of filaments often measured 
before flares /CMEs (e.g., Liu et al. 2012b). We speculate that 
this rising motion signifies the buildup process of the FR driven 
by the tracked converging flows in this AR (Figure 2(c); van 
Ballegooijen & Martens 1989). Nevertheless, the center part 
of the filament is mainly observed to become thickened (see 
Figures 4(b) and (d)) without obvious motion. It is reasonable 
that in the present case, the filament could be supported by the 
magnetic dips below the FR axis (Gilbert et al. 2001), which 
seem to only rise up to ~ 3 Mm (Figure 4(j)). This displace- 
ment is hard to detect considering the projection effect and 
image resolution. 



Figure 5. Temporal evolution of the FR height (a), FR speed (b), maximum a 
near the FR axis (c), total free magnetic energy £f ree (d), and £f ree as a function 
of height (e). The vertical dotted lines indicate the peak times of the HXR 
emission of the X2.1 and X1.8 flares. 

(A color version of this figure is available in the online journal.) 
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1 1/09/06 22:34:22-22: 10:22 1 1/09/07 22:46:22-22:22:22 



Figure 6. Changes of the transverse field in the vertical slice (Figures 4(e)-(l)) from right before to right after the 201 1 September 6 X2.1 (left) and September 7 X1.8 
(right) flares, both showing an obvious leftward (southward) turning of magnetic field predominantly below about 10 Mm. The portion of the slice in gray corresponds 
to the northern end part of the slit drawn in the lower panels of Figure 2, which extends into the outer penumbral region. 


Third, the FR appears to be attached to the photosphere 
before the X2.1 flare; in contrast, it has grown into a semi- 
circular structure with its main body apparently elevated off the 
surface (see Figures 4(e) and (j)). Most interestingly, with the 
release of the nonpotentiality as reflected by a sudden decrease 
of the force-free parameter a (Figure 5(c)), the FR undergoes 
an apparent downward collapse (inward contraction) toward 
the surface after both flares. This sharp evolution represents a 
major restructuring of the coronal field with several noticeable 
features. (1) The spiral of the helical core field lines survives 
the X2.1 flare but vanishes in 1 hr after the XI. 8 flare (see 
Figures 4(g) and (1)). This alludes to the possibility that the 
former event is associated with a partial expulsion of the FR 
(Gibson & Fan 2006; Jiang et al. 2014), while the latter event 
is associated with a fuller FR eruption. (2) The collapse speed 
of the FR is about 1.5 km s -1 in the X2.1 flare and doubles 
to 3 kms -1 in the XI. 8 flare (Figure 5(b)). Thus a more abrupt 
implosion process is very likely during the latter event, which 
is consistent with a larger percentage increase (decrease) of 
the horizontal field (inclination angle) on the photosphere. 
This is the first time that the speed of implosions of the core 
field is derived based on NLFFF extrapolation results. As a 
comparison, Jing et al. (2014) used NLFFF models to find an 
inward contraction of arcade fields overlying a flaring region, 
which had a similar speed (~ 1 .8 km s _1 ) and was also attributed 
to coronal implosions. Here we have to recognize that data 
points during the flares should be interpreted with caution, as 
flare emissions have potential impacts on the field measurement. 
(3) The FR drops in the XI. 8 and X2.1 flares for a distance 
of about 4.2 and 2.1 Mm, respectively. The ratio of 2 in the 
vertical direction is consistent with that (a ratio of 1.6) of 
the decrease of the flux-weighted centroid separation between 
opposite polarities on the photosphere (Section 3.1). Due to the 
vastly different Alfven speeds, however, the time for the coronal 
field to relax to equilibrium takes a few 10 minutes, while that 
for the photospheric field takes 1-2 hr. (4) Our attention is 
also drawn to the change of ambient fields adjacent to the FR. 


In Figure 6, we plot the differences of transverse field in the 
vertical slice between the immediate pre- and postflare states. It 
shows that almost all the fields below about 10 Mm have a clear 
leftward (southward) turning (also see the animation). This can 
be understood since the preflare FR has an inclined orientation 
toward the north due to the pressure imbalance, which might be 
alleviated after the flare energy release. The portion of the slice 
in gray in Figure 6 corresponds to the outer penumbral field 
region, and a more vertical field there following the turning 
nicely agrees with the observed primary decrease of horizontal 
field and decay of penumbra in the Rid and R2d areas (Figure 2, 
lower panels). We conclude that all the above self-consistent 
results make a first comprehensive portray of the 3D coronal 
implosion picture in the low solar atmosphere. 

3.3. Flare Energy Release 

It is also worthwhile to calculate the energy content of the 
homologous X2.1 and X1.8 flares. Solar eruptions are believed 
to be powered by magnetic reconnection which releases energy 
stored in the magnetic field. The upper limit of the available 
energy for energizing flares/CMEs is the free magnetic energy 
£f ree , which can be computed from 



where B N and B p are the modeled NLFFF and potential field, 
respectively, based on the same observed boundary field. For 
all the time instances at which we construct NLFFF, the 
potential field is also derived with a Fourier transform method 
(Alissandrakis 1981) using the vertical component of the HMI 
data. In Figures 5(d) and (e), we present the temporal evolution 
of the volume integrated £f ree and the height profile of £f ree . 
Taking the difference of values averaged over 2 hr before and 
after the flares (denoted by the horizontal red dotted lines), the 
released free magnetic energy A£f ree in the X2.1 (XL 8) flare 
is estimated to be 3.5 x 10 31 (4.6 xlO 31 ) erg. The swift drop 
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Figure 7. (a) Fermi HXR light curves of the 201 1 September 7X1.8 flare. The dotted (dashed) line corresponds to 22:35:52 UT (22:37:22 UT). Spectral fits in the two 
4 s intervals I and II (gray bands) around the HXR peak are shown in panels (b) and (c), using a thermal (green) and a nonthermal (blue) functions. Denoted fitting 
parameters include 8, the electron power-law index, and E c , the low energy cutoff of the nonthermal electron spectrum. The 30-50 keV (between the dotted lines) 
range is excluded in the fitting; see Section 3.3 for details. 

(A color version of this figure is available in the online journal.) 


of £f ree is well cotemporal with the peak times (the vertical 
dotted lines) of the flare HXR emissions (see the next section). 
Before the XI. 8 flare, £f ree accumulates mainly between 3 and 
7 Mm, and decreases at nearly all heights afterward. A similar 
behavior is found before and after the X2.1 flare (also see Jiang 
et al. 2014). 

Another form of energy that we are interested in is the 
energy imparted to nonthermal electrons £nonth- Feng et al. 
(2013) derived a total £nonth of 7.9 xlO 30 erg for the X2.1 
flare using RHESSI data. In Figure 7(a), we plot the lightcurves 
of the XI. 8 flare with the detector nl of Fermi /G BM. This 
detector was the least sunward-facing detector among the four 
Fermi detectors that observed this flare, and it was chosen to 
reduce the pulse pileup effect. Fitting of X-ray spectra was 
conducted for 44 intervals throughout the flare (each with an 
integration time of ~4 s) and in an energy range of ~ 8-300 keV 


(consisting of about 80 energy bins). The 30-50 keV range was 
excluded due to a currently inappropriate handling of the iodine 
K-edge in the response matrix of Fermi ' s Nai detectors. The 
thermal plasma and high-energy emissions are modeled with 
an isothermal (the function vth in SSW) and a nonthermal 
(the function thick 2 in SSW) component, respectively. The 
latter is assumed to be produced by bremsstrahlung emission 
from electrons in a power-law energy distribution (e.g., Emslie 
et al. 2012). Following Ireland et al. (2013), we fix both the 
break energy and the high-energy cutoff to 32 MeV and the 
spectral index above the break energy to 6 whenever appropriate. 
Sample fittings at the HXR peaks are shown in Figures 7(b) and 
(c). Summing the derived nonthermal power over all the fitted 
intervals, the total £ non th of the XI. 8 flare is about 1.6 xlO 31 erg. 

We caution that A£f ree could be up to two times greater than 
the obtained values, mainly due to the preprocessing applied 
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(a) Ale 11/09/07 22:35:52-22:35:07 (b) Ale 11/09/07 22:36:37-22:35:52(c) Ale 11/09/07 22:36:37-22:35:52 (d) AV 11/09/07 22:36:37-22:35:52 



(e) AIA 94 A 11/09/07 22:39:29 (f) TiO 11/09/07 18:15:53 (g) Spectrum of Seismic Emission (h) Egression Power 

1 1 1 1 1 ' 1 1 1 1 ' 1 1 1 ■ 1 i 



Figure 8. Difference images of HMI continuum intensity A I c (a)-(c) and the Doppler velocity AV (d), showing four flare impacts S1-S4 in the 2011 September 7 
X1.8 flare. Selected representative NLFFF lines are overplotted on panel (c). The impact sources are also superimposed as yellow contours on an AIA 94 A image (e) 
and a BBSO/NST TiO image (f), where the blue and red contours depict the identified two seismic sources (at 6 mHz) at the sites of SI and S2. The dotted lines are 
the main segment of the flaring PIL. We show in panel (g) the seismic spectrum of S 1 and S2, and in panel (h) the time profiles of their egression power at 6 mHz. 

(A color version of this figure is available in the online journal.) 


to the photospheric vector data (e.g., Metcalf et al. 2008; Feng 
et al. 2013). Also, E noni ^ should be regarded as the lower limit as 
there exists an uncertainty in determining the low energy cutoff; 
moreover, Enonth of these two flares is derived using two different 
instruments, i.e., RHESSI and Fermi , the systematic difference 
of which is not yet clear. Therefore, the present results can 
only suggest that there seems to be no significant difference in 
the flare energy release between the X2.1 and XI. 8 flares. We 
further note that Fermi /G BM data also show a clear increase in 
count rates above ~300 keV in both flares. A detailed analysis 
of y-ray emission is not intended in this paper. 

3.4. Helios eismic Response to Flares 

To study the seismicity of our homologous flares, we first 
map out flare impacts on the photosphere based on the run- 
ning difference frames of HMI continuum intensity images A I c 
and Dopplergrams AV (e.g., Kosovichev 2011). We then detect 
acoustic emission sources by applying the helioseismic hologra- 
phy technique. Within the scope of this paper, our main concerns 
are the location and relative strength of sunquake sources, and 
their spatiotemporal relation to flare emissions and magnetic 
field structure. 

We begin with the 2011 September 7 XI. 8 flare, which 
showed a seismic activity (Zharkov et al. 2013). As a reminder, 
for this event only HXR time profiles and spectra but no images 
were available from Fermi /G BM. At the initial HXR spike in 
100-300 keV (22:35:52 UT; dotted line in Figure 7(a)), two 
impact sources S 1 and S3 in WL appear near the two ends of the 
flaring PIL in opposite magnetic polarities (Figure 8(a)). Then 
around the HXR peak (22:36:37 UT), SI and S3 become further 
away from the PIL and two additional impacts S2 and S4 are 
formed (Figure 8(b)). According to the NLFFF extrapolation 


(Figure 8(c)), SI is connected to S2 and S3 to S4, which could 
also be confirmed by the 94 A image (Figure 8(e)), with SI and 
S3 lying close to the two footpoints of the erupting FR (see 
Figures 4(c) and (d)). These impacts are also characterized by a 
strong increase of the Doppler velocity up to 1.4 km s -1 (e.g., 
see Figure 8(d)), indicating a downward flow in the photosphere. 
Subsequently, the WL emissions show some separation motion 
away from the PIL as described in the standard flare model, but 
they rapidly become diffused and weakened after the HXR peak 
from -22:38 UT. 

Most intriguingly, despite pronounced WL emission at all 
four locations of S1-S4, it is resolved by the helioseismic 
holography that significant seismic transients are only emitted 
by SI and S2 (also see Zharkov et al. 2013), and that the 
seismic sources are cospatial with S1/S2 at the time around 
the HXR peak (Figure 8(b)). We clearly demonstrate this 
spatial correlation and the source location in the photosphere 
by overplotting the impacts S1-S4 at the HXR peak (yellow) 
and the identified seismic sources (blue and red) as contours on 
a corresponding AIA 94 A image (Figure 8(e)) and an earlier 
BBSO/NST TiO image (Figure 8(f)). To help understand the 
exciting mechanism of the sunquake, we also plot the seismic 
spectra and the time profiles of seismic power of S1/S2 in 
Figures 8(g) and (h), and compare some additional physical 
properties of SI and S2 in Figure 9 and Table 1. Many aspects 
of these results deserve detailed discussions. 

First, several pieces of evidence point to the thick-target 
heating by electrons as a cause of this sunquake event. These 
include that (1) the seismic spectra of SI and S2 have a 
power-law form and display a similarity over a large frequency 
band of ~3-10 mHz (Figure 8(g)). This suggests a common 
source of accelerated electrons from a coronal reconnection 
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Start Time (07-Sep-ll 21:00:00) 


Figure 9. Time profiles of X-rays (a) and some properties of the impact sources SI (blue) and S2 (red) measured at their locations at 22:36:37 UT (Figure 8(b)), 
including the WL intensity contrast (b), red shift (c), LOS field (d), and the vertical Lorentz-force change derived using the vector field data (e) (a negative value 
means a downward force change). The green lines in panel (d) indicate the mean value of LOS magnetic field in the pre- and postflare states. The plotted time period 
of panels (a)-(c) is denoted as the gray band in panel (e). 

(A color version of this figure is available in the online journal.) 


site. These electrons may precipitate into the lower atmosphere 
along magnetic loops linking SI and S2, which are indeed 
both modeled and observed (Figures 8(c) and (e)). (2) The 
seismic sources SI and S2 reach their maximum power almost 
simultaneously at 22:37:22 UT (Figure 8(h)), cotemporal with 
the peak of HXRs (dashed line in Figure 7(a)). This simultaneity 
also implies a direct magnetic connection between S 1 and S2. (3) 


The seismic energy of S 1 and S2 are estimated to be 4.55 x 10 26 
and 1.53 x 10 27 erg, respectively, for frequencies above 3 mHz. 
A stronger seismic emission at S2 could be consistent with the 
fact that the photospheric vertical field at S2 is much weaker 
than that at SI (Table 1), and hence electron precipitation along 
S1-S2 may be stronger at S2 due to the magnetic mirroring 
effect (e.g., Kundu et al. 1995). 
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Table 1 

Physical Properties of Flare Impacts in the 2011 September 7 XL 8 Flare 


Photospheric 

Impacts 

Sunquake 

Source? 

Acoustic Energy 
(10 27 erg) 

Peak of 

A Ic/Ic (%) 

Work by 
SF r a (10 27 erg) 

\B r \ 

(G) 

Horizontal Flows b 
(km s -1 ) 

SI 

Yes 

0.46 

13 

0.36 

880 

0.17 

S2 

Yes 

1.5 

9.4 

0.54 

330 

0.061 

S3 

No 


24 

3.4 

1300 

0.15 

S4 

No 


18 

0.32 

620 

0.12 


Notes. 

a Derived using HMI vector magneto grams. 
b Preflare values. 


Second, the WL emissions at S1/S2 peak at 22:37:22 UT 
(see Figure 9(b)). Thus, they are not only spatially but also 
temporally correlated with the seismic transients. However, 
the impact S2 causes a weaker enhancement of WL intensity 
contrast (9.4%) than SI (13%). Here the image contrast is 
defined as A I c /I c = (/ — h)/h, where I is the intensity of 
the feature of interest and 7o has a value of I in the preflare 
state (at 22:31 UT; e.g., Xu et al. 2006; Kosovichev 201 1). This 
casts doubts on the back-warming of the photosphere, which 
is a secondary effect associated with the WL emission, as the 
mechanism for the sunquake sources SI and S2. We conjecture 
that the initial atmospheric conditions at S1/S2 may allow a 
direct impingement of electrons to cause seismic emissions. 

Third, we also estimate the work Wlf done by the Lorentz- 
force change 8F r on the interior. Following the derivation of 
Alvarado-Gomez et al. (2012), Wlf ~ (l/4)8F r AVAt, where 
AV is the red shift transient seen in Dopplergrams and At is the 
duration of AV . We use two kinds of HMI magnetic field data for 
the evaluation of 8 F r . (1) The LOS magnetograms have the same 
45 s cadence with WL intensity images and Dopplergrams, and 
show step-like changes of LOS field at impacts S 1 and S2 closely 
associated with the flare (Figure 9(d)). As in Alvarado-Gomez 
et al. (2012), we ignore the large excursions around the flare 
impulsive phase that may be affected by flare emissions; instead, 
we average the LOS field over a 8 minute interval both before 
(around 22:30:37 UT) and after (around 22:45: 14 UT) the flare 
(denoted by green lines in Figure 9(d)), and take their difference 
as 8B\ os . For SI in the positive field region, 8B\ os ~ —140 
G indicating a decrease of flux, while for S2 in the negative 
field, 8B los ~ — 62 G indicating an increase. Considering the 
disk location of this AR (N14°W28°) and the orientation of 
the loop S1-S2, the magnetic flux at the SI (S2) region is 
diskward (limbward). Thus the above LOS field changes are 
both suggestive of a collapse of the field lines connecting S1-S2 
(Wang & Liu 2010), which corresponds to a downward Lorentz- 
force change. We derive that the downward force change is 
8F r ^ (1/87T) / dA\8Bi os \ ~ 3.4 x 10 20 dyne at SI and 
1.8 x 10 20 dyne at S2, respectively. With the result of peak AV 
around 22:36:37 UT shown in Figure 9(c) and taking At ~ 90 
s, we have Wlf ~ 2.5 x 10 26 erg at SI and 1.9 x 10 26 erg at S2. 
However, these values would be an underestimation because 
they are only based on the LOS field measurement. (2) We 
also calculate 8 F r using vector magnetograms with a 12 minute 
cadence. The result shown in Figure 9(e) has peaks of 8 F r around 
22:46:22 UT, with which we obtain Wlf ~ 3.6 x 10 26 erg at SI 
and 5.4 x 10 26 erg at S2. We caution that in deriving Wlf above, 
8 F r and AV measured at different times have to be used due to 
the limitation of observations. Nevertheless we do not expect 
this approximation to alter the result significantly, assuming 
that flare-related photospheric magnetic field change tends to be 


stepwise (e.g., Sudol & Harvey 2005). Since for both impacts 
SI and S2, Wlf represents a considerable fraction of the seismic 
energy, the Lorentz-force change may play an important role in 
producing this sunquake. It also needs to be pointed out that 
although 8 F r at S1/S2 amounts only to a few percent of that 
exerted on the region R2 around the PIL (see Figure 2(f)), the 
photosphere undergoes a downward perturbation (as indicated 
by red shifts) predominantly at S1-S4 (Figure 8(d)). 

Fourth, it is obvious that SI is at the poorly developed 
penumbra of the northern sunspot, and S2 is located just at 
the edge of a penumbra segment of the largest spot of the 
AR (Figure 8(f)). In comparison, S3 lies in the umbral region 
while S4 is in the inner penumbral area. In particular, S3 has 
a much stronger vertical field than the other three impacts 
(see Table 1). As suggested by previous studies (e.g., Donea 
2011; Kosovichev 2011), wave motions could be restricted by 
strong magnetic field, which may at least partially explain the 
absence of seismic emission at S3. In fact, seismic sources 
are known to predominantly appear in sunspot penumbrae with 
highly inclined magnetic fields (Donea 2011; Lindsey et al. 
2014). It is however not understood why S4 does not allow the 
development of a seismic source; perhaps, this is related to the 
initial conditions of the low atmosphere. We further notice that 
the region of the impact S2, which induces the strongest seismic 
emission in this event, has significantly weakest vertical field 
and photospheric flow right before the flare (see Figure 2(g) 
and Table 1). We surmise that a weak surface flow field might 
also provide a favorable (but not a necessary) condition for 
the sunquake generation, which needs further observational 
confirmation and interpretation. 

Fifth, it is worthwhile to compare this 2011 September 6 
X1.8 flare with the 2011 February 15 X2.2 flare, both of which 
have multiple sunquake sources and are accompanied by a FR 
eruption. Zharkov et al. (2011) found that in the X2.2 flare, 
two seismic sources are located in the penumbral field at the 
two footpoints of an erupting FR, and that the acoustic emission 
precedes the HXR peak by several minutes and is displaced from 
the strongest HXR source. In the present XL 8 flare, we also see 
that the seismic sources SI and S2 are in the penumbral field; 
however, a noticeable difference is that although the weaker 
SI is cospatial with the northern footpoint, the stronger S2 is 
clearly located away from the southern footpoint of the FR (see 
Figures 4(c)-(d) and 8(f)). The afore-described evidence favors 
the thick-target heating at the HXR peak time caused by the 
direct precipitation of electrons along the loop S1-S2 as the 
sunquake driver in the XL 8 flare. 

After repeating the analysis for the 2011 September 6 X2.1 
flare, we find two bright WL sources (PI and P2) at the HXR 
peak time (see the red contours in Figure 2(c)). However, 
within the detection limit of HMI, no clear seismic sources 


12 


The Astrophysical Journal, 795:128 (14pp), 2014 November 10 


Liu et al. 


Table 2 

Physical Properties of 3D Magnetic Field Restructuring 


Events 

AE free a 

-£nonth a 

Speed (km s 1 2 ) of Implosions 

Distance (Mm) of Implosions 

A B h 

Sunquake 


(10 31 erg) 

(10 31 erg) 

On Surface Vertical Direction 

On Surface 

Vertical Direction 


Association 

2011 Sep 6X2.1 

3.5 

0.79 b 

0.26 1.5 

0.85 

2.1 

26% 

No 

2011 Sep 7X1.8 

4.6 

1.6 

0.17 3.0 

1.4 

4.2 

38% 

Yes 


Notes. 

a Both ALfree and £ n onth are lower limits. 
b Feng et al. (2013). 


or expanding wave ripples can be identified. This is intriguing 
as the X2.1 flare is homologous to the XI. 8 flare, showing 
similar eruption characteristics and energy release. A closer 
look reveals that PI and P2 are mainly around the umbral region 
(see Figures 1(a) and 2(c)), with a vertical field (as well as 
flow field) strength close to that of S3. We thus suggest that 
compared to the XI. 8 flare, the lack of a clear seismic signature 
in the X2.1 flare could be due to the strong vertical field at 
the location of the flare impacts, which may be unfavorable 
for the downward energy transfer, and hence the generation of 
seismic waves, when compared to the inclined magnetic fields 
of penumbral regions (Lindsey et al. 2014). Additionally, a less 
well-formed FR and a slower coronal implosion in the X2.1 
flare may signify a less amount of momentum deposited into 
the surface and below. 

4. SUMMARY AND DISCUSSION 

We have presented a detailed study of the homologous 2011 
September 6 X2.1 and September 7 XL 8 flares concentrating 
on their magnetic field evolution, energy release, and the helio- 
seismic response. By synthesizing the results of photospheric 
field change indicated by observations and coronal magnetic 
field variation suggested by the NLFFF modeling, we have con- 
structed a comprehensive picture of the flare-related 3D mag- 
netic restructuring which turns out to be consistent with the 
coronal implosion scenario in the low solar atmosphere. In the 
aspect of seismicity, we analyzed different flare impacts and 
the associated seismic sources in the X1.8 flare, which provides 
clues to the absence of sunquake during the X2. 1 flare. The main 
results are listed in Tables 1 and 2 and summarized as follows. 

1. In the photosphere, (1) the flux- weighted centroid separa- 
tion of the opposite magnetic polarities shrinks dramatically 
after both flares. We propose that this can be a surface signa- 
ture of coronal implosions. While the approaching motion 
has a somewhat higher speed in the X2.1 flare, the total dis- 
placement in the XL 8 flare is about a factor of two larger. 
(2) The horizontal magnetic field rapidly increases around 
the flaring PIL, but exhibits an asymmetric decrease in the 
surrounding penumbral regions with a much larger magni- 
tude in the northern area. The WL penumbral structure also 
exhibits corresponding changes. Both the mean horizontal 
field strength and inclination angle have larger percentage 
changes associated with the X1.8 flare. 

2. In the corona, (1) the twisted FR lying along the PIL 
collapses toward the surface after the X2.1 flare, then rises 
gradually (up to 0.1 km s -1 ) in 1 day reaching a higher 
altitude, and collapses again after the X1.8 flare. Both the 
distance and speed of the falling motion in the XL 8 flare 
are twice as large as those in the X2.1 flare, indicating a 
more violent implosion and agreeing with a more significant 


change of the photospheric magnetic field. (2) The FR 
is attached to the photosphere before the X2.1 flare; in 
contrast, it is already elevated off the surface before the 
XL 8 flare. A fuller FR eruption could occur in the latter 
event, as little signatures of FRs remain afterward. (3) The 
FR is not symmetric in its central vertical cross section 
but leans northward at 66° relative to the surface. Together 
with the ambient fields, they turn southward rapidly after 
both flares, echoing the observed asymmetric variations 
on the surface. 

3. The released free magnetic energy and nonthermal elec- 
tron energy are evaluated using the coronal field ex- 
trapolation model and HXR observations, respectively. 
Due to uncertainties involved in the computation, it can 
only be concluded that these two homologous flares are 
energetically similar. 

4. Among the four impact sources S1-S4 created by the XL 8 
flare on the photosphere, only SI and S2 spawn seismic 
emissions. Observational and modeling evidence, including 
similar seismic spectra, simultaneous excitation, and foot- 
point asymmetry, all favors the thick-target heating along 
loops S1-S2 over back- warming as the sunquake mech- 
anism. However, the possibility is not excluded that the 
downward perturbation by the Lorentz-force change may 
also contribute. Intriguingly, no sunquake signatures are 
detected in the homologous X2.1 flare with similar erup- 
tion characteristics. A further comparison of magnetic field 
property among S1-S4 and between the two homologous 
flares concur with previous studies that seismic emission is 
more likely to be produced when the flare impact occurs to 
a penumbral region. 

Importantly, our observational and model results portray a 
coherent picture of implosions in the low corona, in which 
the central field collapses toward the photosphere while the 
peripheral field turns to a more vertical configuration, as de- 
picted by Liu et al. (2005). Moreover, the implosion process 
appears to be more abrupt when associated with a fuller FR 
eruption. Compared to the downward Lorentz-force change re- 
sulting from the implosion, the generation of seismic emission 
sources tends to be more closely related to the precipitation of 
the flare-accelerated electrons. Nonetheless, our study demon- 
strates that coronal magnetic field structure and the photospheric 
magnetic field property of the flare impact locations are help- 
ful in understanding the triggering mechanism of sunquakes. 
Extended studies on the spatiotemporal evolution involving 3D 
magnetic field restructuring, particle acceleration, and helio- 
seismic response hold promise of shedding further light on their 
causal relationship. 
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